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ABSTRACT 

Traditionally, secular evolution is defined as evolution of systems where the internal 
growth of structure and instabilities dominates the growth via external drivers (e.g. accre- 
tion and mergers). Most study has focused on "isolated" galaxies, where seed asymmetries 
may represent realistic cosmological substructure, but subsequent evolution ignores galaxy 
growth and interactions. Large-scale modes in the disk then grow on a timescale of order a 
disk rotation period (~ 0.1 — 1 Gyr). If, however, galaxies evolve cosmologically on a shorter 
timescale, then it may not be appropriate to consider them "isolated." We outline simple scal- 
ings to ask whether, under realistic conditions, the timescale for secular evolution is shorter 
than the timescale for cosmological accretion and mergers. We show that this is the case in a 
relatively narrow, but important range of perturbation amplitudes corresponding to substruc- 
ture or mode/bar fractional amplitudes 5 ~ 0.01 — 0.1, the range of most interest for observed 
strong bars and most pseudobulges. At smaller amplitudes S <C 0.1, systems are not isolated: 
typical disks will grow by accretion at a comparable level over even a single dynamical time. 
At larger amplitudes 6 ^> 0.1, the evolution is no longer secular; the direct gravitational evo- 
lution of the seed substructure swamps the internal disk response. We derive criteria for when 
disks can be well-approximated as "isolated" as a function of mass, redshift, and disk stabil- 
ity. The relevant parameter space shrinks at higher mass, higher disk stability, and higher-z as 
accretion rates increase. The cosmological rate of galaxy evolution also defines a maximum 
bar/mode lifetime of practical interest, of ~ 0. If Hubble (z)- Longer-lived modes will encounter 
cosmological effects and will de-couple from their drivers (if they are driven). 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — galaxies: spiral 
— cosmology: theory 



1 INTRODUCTION 

Isolated disk galaxies are prone to a number of important in- 
stabilities that play a major role in shaping observed late-type 
disk and bulge populations, with the most well-known and well- 
studied being the traditional bar and spiral instabilities. Both 
bars a nd spiral structure are ubiquitous in the local disk popu- 
lation dMarinova & Jogeel 120071 ; iMenendez-Delmestre et alj|2007l : 
iBarazza et al. 2008), and their abundance appears comp arable at 
higher redshifts dsheth et al.l20oil2008l:|jogee et al.l2004l) . By am- 
plifying small perturbations into coherent, long-lived, large-scale 
non-axisymmetric modes, these structures enable disks to evolve 
significantly - re-distributing material in angular momentum and 
phase space - in a few orbital periods. As a consequence, observa- 
tions and simulations indicate that these structures are important in 
shaping the cosmological evolution of disk sizes, scale heights, and 
the abundance, structural properties, and mass fraction in "pseu- 



* E-mail:phopkins @ astro.berkeley.edu 



dobulges" (disk-like bulges that result from angular momentum 
exchange in these modes), a population inc reasingly prominent in 
low-mass and later-type disk galaxies (e.g. Debattist a et al.ll2004l : 
iKormendv & Kennicuttll2004iweinzirl et alj|2009l) ~ 



Traditionally, the growth and evolution of these global modes 
is referred to as "secular" evolution: by definition, evolution that is 
slow relative to the local dynamical time. This contrasts with vi- 
olent relaxation - seen in e.g. galaxy-galaxy major mergers - in 
which the potential fluctuates on short timescales, and local insta- 
bilities, involving e.g. clumping, star formation, and formation of 
bars on small scales (sub-kpc). 

As a consequence, the secular evolutionary channel has, for 
the most part, been studied in the context of isolated galax- 
ies. Given an isolated, self-gravitating stellar (or stellar+gas) disk 
that meets certain instability criteria, small non-zero amplitude 
in the large-scale modes that identify morphological bar and spi- 
ral patterns (characteristic wavelength of order the disk length) 
will grow exponentially on a timescale of a few orb ital periods 
(see e.g. the discussion in iBinnev & Tremainel Il987t) . The evo- 
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lution and dynamics of these modes have been well-studied in 
idealized cases of isolated disks with properties simil ar to the 
Milky Way, but by design bar or spiral wave-unstable |Schwarj 
Athanassoula et al 



1981; 



Comb es et al 
19931: 



19831; lPfennigeri[19 84; Wein berg! 1 19851; 
1990; Hernquist & Weinberg 1992a; Friedli & Benz 



Athanassoula 
2002b; 



I2002al; lAthanassoula & Misiriotij |2002|; 
Weinberg & Katz 2007b; Kaufmann et al. 
' 20001; iMaver & Wadslevl |2004 



Athanassoula 

2007; Patsis & Athanassoula _ . . 

Berentzen et al.ll2003L |2004| ; iFovle et alj|2008» . In particular this 



informed study of the role of secular evolution in shaping galaxy 
sizes, dynamics, and morphology. 

However, in ACDM cosmologies, structure grows via contin- 
uous accretion and mergers. Although major mergers are rare, both 
theoretical calculations and observations suggest that minor merg- 
ers are ubiquitous, and accretion of new cold gas is rapid in low- 
mass galaxies (Woods et al. 2006; Mail er et al.ll2006l ; lBarton et al.1 
120071 ; Iwoods & Gellej|2007l ; IStewart et alJl2008bl)~Togefher with 
the typical substructu re present in ACDM halos dTavlor & Babull 
l2004l ; lGao et alj2004h . this suggests the concern that there may not 
be in practice such a thing as an "isolated" galaxy at the level of 
interest. 

More recent studies of secular evolution have therefore fo- 
cused on more realistic scenarios, exploiting merger histories from 
cos mological simulations in semi-idealized studies of single galax- 
ies (Bournaud & Combesll2002l;lBenson et"ai]|2004 Gauthier et al. 



20061; Berentzen & Shlosmanl 



200' 



Curir et alj 120071 ; iKazantzidis et alj 120081 ; iRomano-Diaz et al.l 



Kaufmann et al 



20071; 



20081) . These simulations again reveal bars and spiral structure 
to be prominent - arguably more so than in isolated simulations 
- but it is less clear whether their formation and evoluton can 
be attributed to the same secular processes at work in isolated 
systems, or whether they are driven systems owing to substructure 
and accretion in the galaxy disk and halo. 

The important question for models is: can any galaxy in a re- 
alistic cosmological context still be approximated as "isolated" for 
certain purposes? If so, in what regimes as a function of redshift, 
galaxy mass, and internal properties is this applicable? What are the 
corresponding implications for interpretation of bar fractions and 
lifetimes? And ultimately, what does this imply for the importance 
of isolated secular evolution in driving the evolution of galaxies and 
formation of bulges? 

In this paper, we attempt to address these questions by means 
of a simple comparison of cosmological accretion rates and charac- 
teristic timescales for secular evolution. This approach allows us to 
identify the regimes where galaxies can be safely considered "iso- 
lated" versus where cosmological effects may not be negligible. 
We show that there is an interesting regime of secular modes with 
fractional mass/amplitude ~ 0.1 where the secular growth mode 
dominates and the isolated galaxy approximation is good (§[2j». We 
show how this scales with galaxy mass, redshift, and disk stability 
properties (§[3](, and identify some basic consequences for the life- 
times of large-scale modes in disks (§[4]l. Our goal is not a definitive 
description of secular evolution, but rather to provide a set of sim- 
ple initial constraints to provide context for more detailed studies 
of the interesting parameter space. 

Throughout, we adopt an Qm = 0.3, J1a = 0.7, h — 0.7 cos- 
mology, but our conclusions are not sensitive to th e choice within 
the r ange allowed by present observations (e.g. iKomatsu et al.1 
12009b . 



2 SECULAR EVOLUTION VERSUS COSMOLOGICAL 
EVOLUTION 

Consider an "initial" equilibrium, axisymmetric disk+halo system 
at time t — 0. In this limit the system will not evolve any non- 
axisymmtric modes. Therefore, introduce a non-axisymmetric per- 
turbation to the disk potential of amplitude 



So = ^- 



(1) 



We are specifically interested in global models, so 4> ~ GM/R is 
the potential of interest (where M is the disk+enclosed halo mass 
and R is a characteristic effective radius/scale length). The precise 
meaning of the perturbation 8cj> differs depending on the mode(s) of 
interest and configuration. For example, in idealized N-body simu- 
lations, this typically corresponds to shot noise. However, in realis- 
tic cosmological settings this will correspond to substructure in the 
disk or halo, with 5(j> ~ Gm/r (where m is the substructure mass 
and r its "initial" distance). The relevant numerical prefactor will 
depend on the orbit, phase-space structure, and mode (for example, 
for a bar, the desired quantity is the time-averaged contribution to 
the m = 2 mode at radius ~ R in the co-rotating frame); for our 
purposes, the scaling (not absolute value) of 8 is most important. 

At early times (before saturation), this non-axisymmetric term 
will be amplified internally and grow roughly exponentially: 



S(t) = So expO/fo) 



(2) 



where to is the effective secular timescale, which is typical of order 
a fe w orbital times (again, this is for global modes, not local; see 
e.g. iHollev-Bockelmann et all 120051: [Weinberg & Katzl l2007al lbh. 
This growth time has been the focus of considerable study, and 
is one of the many important results of isolated disk studies. For 
example, for a disk bar in a strongly unstable bulge-free Milky 
Way-like disk, Dubin ski et alj J2008) show Equation l[2j is a good 
approximation to the behavior in simulations, with to = 9>tt / k ~ 
2.83 Pd (where k is the epicyclic frequency, = 2 3/ ^ 2 7rP rf ~ 1 for a con- 
stant circular velocity disk, and Pj = 2nR/V c is the disk circu- 
lar period at its effective radius). iKlypin et alj |2008) find a sim- 
ilar fo ~ 3 — 5 Pj for thin, bulge-free MW-like disks (albeit with a 
much larger to ~ 10 — 30Pd for thick H/R j> 0.5 disks; see also 
IColm et al1l2006l) . iMartinez-Valpuesta et alj j2006l) see timescales 
from ~ 2.5 — 10 Pd, depending on whether the bar growing is an ini- 
tial mode or a secondary (post-buckling) mode. A similar range of 
timescales is found (with co nsiderable galaxy-to-galaxy va riation) 
in live cosmological halos in lBerentzen & Shlosm an (2006). 

For less cosmologically motivated, but more general and an- 
alytica lly tractable disk mass profiles, I Athanassoula & Sell wood! 
( 1986) find typical to ~ 1 .0 — 6.7 Pd for realistic halo mass fractions 
~ 1/4—1/2 (fraction of the t otal mass owing to th e ha lo at < R) 
and sc ale heights H/R ~ 0. 1 . iNaravan et al.l dl987h and lShuet ail 
jl990h obtain to ~ 0.8 — 1 Pd for gas disks with an outer Lindblad 
resonance alR>R e (of interest for global modes here) and no halo. 

More stable systems will evolve more slowly; for the sake of 
generality we define 



to = Afoisk Periods X Pd = 



1 



Xeff 



-ft 



(3) 



where Xets is an effective stability parameter: Xea ~ represents 
typical, cosmologically realistic disks maximally unstable to large- 
scale modes, which will evolve on a single orbital time; and XsB > 1 
systems are stable and experience only oscillations, rather than 
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amplifying modesQ Note that, formally speaking, Xcff < is al- 
lowed. For certain bar configurations, for example, to ~ 0.7 fi; ha s 
been obtained (see e.g. lAdams et al. 1989: E arn & Sellw ood 1995), 
or even, for spiral structur e in the weak winding approximation, 
to ~ 0AP d jToomre|[T98lh . However, those situations all involve 
no halo and an infinitely thin disk, and somewhat different matter 
profiles from what are observed in typical disks. For moderate halo 
contributions or disk thickness, to is unlikely to be smaller than Pj 
by any but a small factor (to ~ 0.7 — 0.8 Pj), a small difference rela- 
tive to the uncertainties in other quantities calculated here. And the 
MW-like examples above illustrate that Xeff ~ 0.5 — 0.75 is proba- 
bly the case of greatest interest for realistic disk plus halo systems, 
even for strongly unstable systems. To be conservative, however, 
we will adopt Xcff = for all numerical estimates, unless explicitly 
otherwise specified. 

However, galaxies are not static, and two things will happen 
that might compete with this internal self-amplification: (1) the sub- 
structure itself can dynamically evolve, driving stronger perturba- 
tions and/or merging; and (2) new mass of magnitude comparable 
to the disk mode can be accreted/merge. If either of these occurs 
on a timescale shorter than to (the effective secular timescale), the 
system should not be considered "isolated" for purposes of secular 
evolution. 

Consider case (1), the dynamical evolution of the substruc- 
ture itself. Given some substructure/perturbation of mass fraction 
5 at some initial radius of interest r, the orbit will decay on a 
timescale of order the dynamical friction time; correspondingly, 
the perturbation 5 oc 5<f> oc r _I will grow on the same timescale|j 
Strictly speaking, dynamical friction does not dominate angular 
momentum loss at small r adii; rather, resonant tidal interactions 
act more efficiently dBarnes & Hernquistlll992h . However, prop- 
erly calibra ted, the dynamical friction t ime is not a bad approxi- 
m ation (e .g. Boylan-Kolchin et al. 2008!). For an isothermal sphere 
or lMestelldl963h (flat rotation curve) disk, this time is simply: 



fdf = 



R/Vc M mc (r) r 
2/3 In A m R 



0.2 
8 InA 



Pa 



(5) 



' Under certain restrictive circumstances, our Xcff h ere is analogous to 
the Toomre Q or X parameter X = k 2 R/ (2nnGTi) or a (renormalized) 
Ostriker-Peebles criterion (proportional to the ratio of rotational kinetic to 
potential energy). For e.g . the bar in a two-dimensional Ku z'min disk ap- 
proximation presented in Athanassoula & Sellwoodl il986l> . we can trans- 
late their Equation 3 to obtain 



Xeff = 0.3+1.1 



' /halo + Aulj, 



1 ) +0.6 ( 



H 
OAR 



1 



(4) 



in physical terms of the disk thickness H/R and halo plus bulge (non-disk) 
mass fraction inside R. The definition in Equation[3]is not, however, meant 
to represent specific instabilities, but to allow for general large-scale disk 
modes with a characteristic growth time/stability criterion. 
2 Strictly speaking, realistic cosmological perturbations grow continu- 
ously, so an "initial" radius is ambiguous. However, there is still some 8<f> 
that scales as described at a given instantaneous r, and this is what ulti- 
mately enters into the equations derived. Also, in practice, such modes - 
where induced by substructure - often appear suddenly (i.e. in a time < Pj 
when r ~ R; this is because at larger radii, the net non-axisymmetric 8</> 
contribution is suppressed by a Poisson ~ N~ 1 / 2 (~ R~ 3 / 2 ) term. In sim- 
ulations, for example, perturbations are typically dominated by a few close 
passages of clumps/substructure where r ~ R (althou gh these may be from 
longer radial orbits; see Velazauez & White 1999; Bournaud & Combes 
l2002l:lGauthier et al j|2006l : iKazantzidis et al jhoOSt iHopkins et alj|2008bl) . 
In any case, since our derivations rely on 5, rather than r explicitly, this is 
not a large source of uncertainty. 



where the equality on the right comes from the definitions of So and 
Since we are considering the magnitude of the perturbation 
relative to the disk, the time here scales with the disk dynamical 
time at fixed So (as opposed to e.g. the Hubble time for halo-halo 
orbital decay at large radii). 

The left panel of Figure[T]compares this timescale to the sec- 
ular evolution timescale to. For representative purposes, we as- 
sume a "maximally unstable" to = Pi (Xcff = 0) MW-like disk 
withP rf = 27r x 5kpc/200kms~' w 160 Myr, and total stellar mass 
= 5 x 10 10 Mq. This is easily generalized; Pj (at the scale ~ R 
of the disk itself) appears to be independent of mass in observed 
disks (e.g. ICourteau et all 120071 . and references therein). We plot 
the results assuming such a disk exists at redshift z = 1, but the 
qualitative scalings are similar at redshift z — 0, and we will show 
the redshift dependence explicitly below. We compare the dynami- 
cal friction time tat; here we show the result s using the full orbital 
param eter-dependent fits from simulations in lBovlan-Kolchin et al.l 
(2008), which allows us to quote the ±la range of fdf from 
the range of orbits observed in c osmological simulations dBensonl 
120051 : iKhochfar & BurkerJ 120061) . Using the simpler formula in 
Equation[5]is similar to the median expected. 

Comparing Figure Q] or Equations [3] & [5] shows that the dy- 
namical evolution of the perturbation is more rapid than the internal 
response for mass ratios larger than 



<5crit,df = 



1 



4 7T /3 InA 



'0.2(1 -Xeff) 



(6) 



This is ultimately an obvious regime; when 8(j>/4>^ 1, direct evolu- 
tion dominates the potential fluctuations. We denote this the "major 
merger regime": in the case where 8 corresponds to some substruc- 
ture, this clearly requires a mass ratio /j, > 0.2 with r^Rj, i.e. close 
passages of major companions. Note though that this does not have 
to be a merger. For example, a sufficiently strong disk fragmenta- 
tion event will be similar. Physically, this is still dynamically dis- 
tinct from secular evolution (from e. g. bars, etc.) - it will "look 
like" a merger inside the disk (see e.g lElmegreen et al]|2008l) . 

Now consider case (2): new growth/perturbations/mergers. 
Note that we are no longer considering the evolution of indi- 
vidual perturbations, but the time between new perturbations of 
the same or greater magnitude. If this is -C to, then the sys- 
tem is not isolated. A lower limit to this is given by the rate of 
baryonic accretion/merging onto the disk (if accreted systems re- 
tain some dark matter, they will represent larger perturbations, 
but there is at least a lower-limit in the mass added in baryons 
to explain the disk mass). Detailed analyses of thes e rates have 
been discussed extensively in the literature (s e e e.g. | Brown et all 
20071; iGuo & White! 1200a IWetzel et al.l 120091 : iGenel et alj|2008l: 



Stewart etall20 08b). Here, we use a simple semi-empirical model 
to define som e of the relevant scalings; for more discussion, see 
Hopki ns et al] d2009al) . A variant of the model, bas ed on subhalo- 
subhalo merger rat es, is also described in detail in lHopkins et al.l 
(2008a). Following Stewart et al. ( 20 08*3), we begin with dark mat- 
ter halo merger trees (here from Fakhouri & Ma 2008). Empiri- 



3 In detail, /3 is a constant that weakly depends on the mass profile and 
velocity is otropy: = 0.42 8 for an isotropic isothermal sphere and = 0.32 
for a thin iMestell 1 19631) disk averaged over random inclinations (used 
in Equa tion [3J. Th e Coulomb logarithm is a pproximately A = 1 + 1 /So 
I Boylan-Kolchin et al. 2008; Jiang et al. 2008). For FigurefT] we use the fit- 
ting functions from IBovlan-Kolchin e t al. (2008), with appropriate eccen- 
tricity and orbital parameter dependence, rather than the simplified Equa- 
tionfS] but the results are similar on average. 
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Figure 1. Characteristic timescales for evolution of perturbations in unstable ~ L» disks (here at z = 1, with Xcff = 0). Left: Timescale for internal disk 
response (secular evolution) to amplify some large-scale mode with amplitude S, compared to the timescale for an individual perturbation to evolve on its own 
(via e.g. dynamical friction). Analogous to major mergers, direct evolution is more rapid than disk response for major perturbations > 0.2. Units are the age 
of the Universe at this redshift. Center: Secular timescale versus timescale for the disk to accrete a new fractional gas mass > <5 or undergo a new merger with 
mass ratio > 6. At sufficiently low amplitudes, accretion is non-negligible over the secular response timescale. Right: All timescales. Disks are effectively both 
isolated and potentially secular evolution-dominated in a regime around 8 ~ 0.1. Raising Xeff will increase the "secular timescale" and decrease this range. 



cal halo occupation models and other observations constrain the 
average galaxy mass per host halo (or subhalo) mass, with little 
scatter - so at a given instant we simply populate the halos with 
galaxi es. Here assigning s t ellar m ass given the fitted M* (Mhaio | z) 
from IConrov & Wechsleil <l2009t) and gas mass given the fits to 
M gas (M* | z) from lStewart et al.l ?2009) (for the observations used in 
the fi t s, see references therein and [ Bell & de Jondl200ll : lErb et al] 
120061: iFontana et alfeOO&IPerez-Gonzalez et al]|2008h . 

The uncertai nties in this modeling methodology will be dis- 
cussed in detail in Hopkins et al. (2009b), but for our purposes they 
are relatively small (factor ~ 2 uncertainty in the merger rate near 
~ L„, owing to a combination of uncertainty in M» (Mhaio) and the 
halo-halo merger rate) at z < 2, because it is primarily the shape 
of the galaxy-halo mass correlation (rather than e.g. its absolute 
normalization) that affects galaxy-galaxy merger rates Q Note, how- 
ever, that the uncertainties grow rapidly at higher redshifts, owing 
to the lack of empirical constraints. Evolving the system forward 
some small increment in time, we can "add up" the mergers (in de- 
tail, we add a dynamical friction "delay" time between each halo- 
halo merge r and subsequent galaxy-ga laxy merger, with the for- 
mulae from Bov lan-Kolchin et alj|2008l) . This gives merger rates; 
but also, knowing the new halo mass (after accretion/growth in this 
time interval), the empirical halo occupation constraints define the 
"expected" galaxy mass for the updated halo mass. We simply as- 
sign whatever galaxy mass growth is needed to match this (not 
already brought in by mergers) to "accretion." Note that this is a 
lower limit to the accretion rate, reflecting net accretion (outflows 
may remove mass, requiring more new gas inflow). 

The middle panel of Figure Q] shows the relevant timescale 
for both mergers (median time At between mergers with baryonic 
mass ratio fi = Mbar^/Mxir. i > <5o) and accretion (At for the disk 
to grow via accretion by a mass fraction > <5o). Accretion tends 



4 The merger rates from this model as used here can be also ob- 
tained as a function of e.g. galaxy mass, mass ratio, and red- 
shift from the "merger rate calculator" script publicly available at 
http : / /www .cfa. harvard. edu/ ~phopkins/ Site/merger calc 



to be the dominant growth channel (relative to e.g. minor merg- 
ers), for all but the most massive galaxies (where gas accretion is 
"quenched"). As a result, the time between new mergers may be 
long, but at sufficiently low So, growth by accretion is more rapid 
than internal disk response. We denote this the "accretion regime." 
Again, the behavior is easily understood: if one is interested in evo- 
lution at the <C 10% level, then galaxies cannot be considered iso- 
lated for even a single dynamical time, as they will grow by more 
than this amount in that time. 

The relevant criterion can be roughly estimated as follows: to 
very crude approximation, fractional galaxy growth rates scale as 
~ a/fHubbie, where a is weakly redshift dependent but non- trivially 
mass-dependent with a ~ 0.2 for a Milky- Way mass halo at z = 1 
(i.e. an assumed galaxy mass 5 x 10 I0 Mq). For such a system, as 
pictured in Figure[T] the galaxy will grow by a fraction > So in the 
time to (secular response time) for perturbation amplitudes below 



Pel 



0.003 



u crit. acc — 



(z = 0) 



(7) 



(1-Xeff) t H (z) 1 -Xeff 

Figure Q] considers the "maximally unstable" (xeff = 0) case, 
such that to = Pj. If the stability parameter is higher (larger to), 
the regime of effective "isolation" will be more restricted. Figure|2] 
illustrates the parameter space as a function of the effective disk 
stability parameter Xets (recall, this is simply defined relative to the 
number or orbits needed to grow the mode of interest). Above some 
critical Xeff (here Xcff ~ 0.75, i.e. Ambits = 4 or f o > 0.5 Gyr for a 
MW-like disk), the secular timescale is always longer than the other 
timescales above. This is simply the statement that disks are not 
"isolated" for timescales > Gyr, especially at high redshift. 



3 DEPENDENCE ON GALAXY MASS AND REDSHIFT 

Figure [3] shows how the regime of secular evolution depends on 
galaxy mass and redshift. First, we consider the same compari- 
son at z = as a function of galaxy mass. Observations indicate 
that Pj is nearly mass-ind ependent at the disk effective radii of in- 
htterest for global models (Bell & de Jong 2001; Shen et al. 2003; 
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Figure 2. Parameter space of regimes in Figure □for the same ~ L* system 
as a function of perturbation amplitude and effective disk stability (speed of 
the growth of secular modes of interest). 



ICourteau et alj|2007h . Given Equation [5] the same is true for dy- 
namical evolution of individual perturbations (at fixed So). How- 
ever, accretion and merger rates scale significantly with mass. At 
low masses, merger rates are low, but accretion rates are high. 
At high masses, accretion rates drop rapidly (consistent with zero 
at M„ai 2> 10 1 Mq), but merger rates increase, leaving almost 
no range of perturbation in which secular processes are relevant 
(right column of Figure [3). Both effects are seen in a variety of 
models and observations ([M ailer et a l.l 120061: iNo eske et al. 2007; 
Guo & Whitell200ot iKitzbichler & Whitdl2008l; IParry et aljbooa 
et alj l2008al : iKeres et al.l 120091 : [Bundv et alj 120091) . The 



Stewart 



mass dependence is important even over a relatively narrow mass 
range - for example, note that our previously assumed Milky-Way 
like mass of 5 x 10 10 Mq (Figure [TJ, being a factor ~ 2 smaller 
than the 10 11 Mq case shown here, has correspondingly more rapid 
accretion rates (between the 10 lc Mq curve and 10 11 Mq curve). 

Again, we emphasize that we are using baryonic mass ratio fi 
here - this is a minimum, as it reflects the most densely bound ma- 
terial that will survive to perturb the galaxy (an individual merger 
may "begin" at larger 5 including dark matter, or smaller S at large 
radii, but orbital decay and stripping will tend to saturate it at 
54>/4> ~ with a rate of new su ch events from mergers as shown; 
see e.g. iKazantzidis et alj|2008l) . Low-mass galaxies are observed 
to be more dark-matter dominated, so if this can be conserved, the 
relevant rates will not decrease as rapidly with stellar mass; how- 
ever, modeling this requires more detailed knowledge of cosmolog- 
ical orbits, stripping, and internal galaxy structure. 

For each mass, Figure [3] shows how the regime of secular evo- 
lution depends on redshift. To lowest order, accretion timescales 
evolve with the Hubble time (fitting directly, accretion rates oc (1 + 
z) 2 ; see lStewart et alj|2008al) . Observations of the baryonic Tully- 
Fisher and size mass relation suggest that Pj (or eq uivalently at 



fixed mass, disk sizes) evo lv es weakly from z = 



e quivalently at 
2 dTruiillo et al.1 



20061: iFlores et alj 120061: iKassin et al] 120071; iToft et al.1 |2007t 



Akiv ama et alj2008l ;l Somerville et al .120081) . Moreover, theoretica 



models that include the well-establis hed dependence of halo con- 
centra tion on redshift (see e.g. lBullock et al. 2001; Wechsl er et al.1 
120021) predict a similar weak scaling dSomerville et al. 2008]). Pa- 
rameterizing as Pd oc (1 + z)~ l3d , these observations constrain /3d = 
0.0 — 0.6. In Figure(3] we conservatively adopt = (i.e. Pd inde- 
pendent of redshift), but we show how the results would change 
if we allowed the maximum observationally inferred evolution, 



f3d = 0.6. It makes a small difference, but does cancel some of the 
redshift evolution in the relevant parameter space. Even in the ex- 
treme case of a simple Pd oc fHubbte scaling dMo et al.l l 1998). some, 
but not all of the evolution is negated (at z < 2, merger and cooling 
rates evolve as oc (1 +z) 2 , l/tnme as oc (1 +z)). 

In Figure[3] the critical amplitude below which the "accretion 
regime" pertains scales roughly as Sah, acc oc (1 + z) L5 ~ 2 ' , while 
5crit,df ~ constant. This is an approximation over the entire range 
z — — 2; in fact at the lowest redshifts (z < 0.2), the falloff in 
<5 C rit, acc is somewhat more rapid (as e.g. the Universe's acceleration 
term becomes important). As a consequence, the range of So over 
which "isolation" is a good approximation decreases with increas- 
ing redshift. 

Figure [4] summarizes the parameter space as a function of 
galaxy mass and stability parameter XeB, at z = 0, z — 1, and 
z — Zform(M g ai). We define Zform(M ga i), the galaxy assembly time, as 
the redshift when each galaxy reaches half its z = mass, accord- 
ing to our simple growth model. To the extent that secular modes 
are considered important in this formation process, this is an inter- 
esting timescale. 

Simulations find that star-forming galaxies accrete most of 
their mass along a couple of dynamically coherent, clumpy fil- 
aments; as s uch they are dynamically important for large-scale 



aments; as sucn they are dynamically important tor large-scale 
disk modes dKeres et alj 12001 120091 ; iDekel & Birnboiml 120061: 
iDekel et alj|2009l) . If, however, accretion were perfectly smooth, 
axisymmetric, and restricted to large radii (without migration of 
new material inwards), then it might be valid to ignore it in study- 
ing secular modes even when accretion rates are large. To represent 
this possibility, Figure [5] re-calculates Figure [4] but ignores accre- 
tion. At low masses, merger rates are sufficiently low that the iso- 
lated regime extends to smaller mass ratios 5 < 0.01. 



4 IMPLICATIONS FOR MODE "LIFETIMES" 

The cosmological evolution of galaxies also has important impli- 
cations for mode "lifetimes." Since frj = 1/(1 — Xes)Pd, there is 
clearly some Xcff at each redshift above which to is larger than any 
of the competing timescales for all So- Modes with larger Xes are 
still formally unstable, but the time/number of orbits to amplify 
the mode becomes sufficiently long that these modes should be 
considered cosmologically dynamical objects. Figure [6] shows this 
maximum Xtff as a function of redshift (for ~ 10 u Mq galaxies 
where this is maximized, as seen in Figure[4]l. At z > 1, this cor- 
responds to modes growing in < a couple Pd\ at z > 2, however, 
even Xcff =< 1 systems (those where modes grow on a timescale 
~ Pd) can be in the "accretion regime," as discussed above. Recall, 
simulations suggest that eve n cold, bulge-free MW -like disks have 
effective Xcff ~ 0.5 — 0.75 dDubinski et alj I2OO8L and references 
therein). This high-z behavior is directly related to observations 
showing that disk orbital periods at high redshifts become com- 
parab l e to the Hubble t i me (see e. g. Flores et al. 2006; Kassin et al 
2007 1 ; lToftetal.ll2007l; Ivan Starkenburg et alj|2008t Ishapiro et~ 
2008). 



At Xeff less than the values above, secular modes can grow 
"in isolation" from some So. Typically, these will grow rapidly and 
saturate at some 5f ~ 1. However, if an isolated mode then sur- 
vives stably at an amplitude 5f for a lifetime much longer than 
the other timescales compared here, then various cosmological ef- 
fects may have important consequences. Fo r example, if a disk ba r 
saturates and survives with some 8j ~ 0.4 toubinsldetaiT 2008), 
in some number of dynamical times the galaxy will grow by this 
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Perturbation Amplitude 5 



Figure 3. Same as the right panel of Figure [T](Xcff = 0), as a function of galaxy mass and redshift. In low-mass (< 10 10 Mq ; left) galaxies, merger rates are 
low, but accretion rates rapid - secular responses at the <C 10% level compete with cosmological disk growth. At intermediate masses ~ L, (~ 10" Mq; 
center) accretion and minor mergers occur with comparable rates. At high masses (> 10' ' Mq ; right) accretion rates are low (cooling is inefficient) but merger 
rates grow rapidly - secular responses at the < 10% level compete with mergers. As a function of redshift, disk dynamical times scale weakly, but merger 
and accretion rates increase, leaving less of parameter space in which disks can be considered "isolated" for the internal response time. Error bars mark the 
range between the internal response time if disk sizes do not evolve with redshift (fit = 0; dashed lines) and if they evolve at the maximum rate constrained 
by observations (J3j = 0.6; lower bar). 



much. Essentially, cosmological growth may "catch up" to the sat- 
urated mode and could effect it. Of course, the mode could continue 
growing with the galaxy, or be robust to these effects; our point is 
that continuing to treat such a mode in isolation may not necessar- 
ily be a good approximation over much longer timescales. More- 
over, if stable modes can survive for a timescale much longer than 
e.g. the relevant dynamical friction times at So, then the presence 
of those modes mode (the duty cycle) will de-couple from that of 
their drivers (if they were initially driven). In e.g. the case of minor 
mergers, this is the statement that new mergers and/or the destruc- 
tion of the original driving satellites will wipe out the "memory" of 
the drivers, while the bar survives. 

Taking the minimum of the non-secular timescales of inter- 
est (e.g. accretion and merger timescales in Figure [3} at whatever 
amplitude 5 maximizes this timescale, gives the maximum relevant 
"isolated" mode lifetime. This is clearly a function of mass; we 
consider here the ~ 10 1 1 Mq (~ L*) case, of greatest interest both 
as a MW-like system and because Figure[4]demonstrates that this is 
where such a timescale (the "isolated" regime) is maximized. Fig- 
ure |7] plots this timescale versus redshift. We show this both for 
the assumption that Pj does not evolve (/3j = 0) and the maximum 
observationally constrained evolution (f3j — 0.6). We compare a 
constant fraction (~ 0.1) of the Hubble time - this appears to be 
a good approximation, on average (there will of course be scatter 
galaxy-to-galaxy in accretion and merger rates, leading to typical 
factor ~ 2 scatter in the relevant timescale here). 



5 DISCUSSION 

Under typical cosmological conditions, global "secular" evolution 
- narrowly defined as evolution by internal amplification of large- 
scale disk modes in effectively isolated galaxies - only occurs in 
a restricted range of parameter space (Figures [U2t . If the pertur- 
bation mode of interest has a fractional amplitude <C 0.1, what we 
call the "accretion regime," then the disk will grow by accretion by 
a comparable amount in even a single dynamical time; the isolated 
approximation is clearly not valid. This threshold is around an am- 
plitude dcm, acc ~ 0.002(1 -Xeff) -1 (1 +z)' - 5 ~ 2 for 1O U M galax- 
ies (slightly lower at z < 0.2) or Sent, acc ~ 0.005 (1 - Xeff) - ' (1 + 
z) 15-2 for 10 10 Mq systems. At the opposite extreme, seed "per- 
turbations" of fractional amplitude So > <5 C i-it, dt ~ 0.2 lead to non- 
secular evolution - the perturbations' own gravitational evolution 
will dominate the internal response (this is obvious in the case of 
e.g. galaxy-galaxy major mergers or massive disk fragmentation 
events, where the evolution of the merger/clumps drives violent re- 
laxation). 

The relevant parameter space depends on galaxy mass (Fig- 
ures \315\. Although halo growth is nearly mass-independen t 
dFakhouri & mH|2008| ; IGuo & W hite 2008; IStewart et alj|20083) . 
galaxy growth histories are not (the function M ga i(Mh a io) is non- 
trivial). At high masses (M gll i > 10 u Mq) galaxy-galaxy merger 
rates are high such that systems are rarely "isolated" over the 
timescales of interest for secular evolution. However, there are few 
disks at these masses, so secular evolution is not expected to be 
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Figure 4. Parameter space of the "isolated" secular regime (Figure[TJ versus 
galaxy mass and Xcff. for a redshift-independent disk period Pj. Top: At the 
formation redshift z = Zf or m> where each galaxy reaches half its z = mass. 
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Figure 5. As Figure |4] (with z = Zform); but neglecting cosmological accre- 
tion (i.e. considering only merger timescales as the competing timescale at 
low-<5). 
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maximum isolated lifetime in Figure [7] Black (red) lines show the median 
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Figure 7. The maximum "isolated" lifetime of large-scale modes (e.g. disk 
bars). This is the largest timescale (marginalizing over <5 and Mg a i) shorter 
than the competing (non-secular) timescales in e.g. Figure [3] Variations 
within populations contribute factor ~ 2 scatter from object-to-object. We 
show a fixed fraction of the Hubble time for comparison. Evolution on 
longer timescales will compete with cosmological effects. The duty cycle 
will de-couple from driving: even if all such modes were driven by e.g. mi- 
nor mergers, there will be no correlation between the presence of the modes 
and the presence of companions. 



a dominant process. At low masses (< 10 10 Mq) merger rates are 
low (in terms of galaxy-galaxy baryonic mass ratios; including dark 
matter, they may remain high) but accretion rates are high; systems 
can be effectively approximated as isolated for only a couple of or- 
bits in the regime of amplitudes 8 ~ 0.03 — 0.2. Moreover, although 
such galaxies are mostly disk (B/T <C 1), they are increasingly 
dark matter-dominated which hel ps stabilize them t o the devel- 
opment of secular modes (see e .g. Persic et al. 1996; Miho s etakl 
1 19971; iBorriello & Saluccill200ll; Bell & de Jondl200lh . Galaxies 
may be "most isolated," and so traditional secular evolution most 
relevant, between these regimes, i.e. in galaxies somewhat below 
~ L» . That this occurs at masses only somewhat below where merg- 
ers become efficient is also interesting; there may be a relatively 
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rapid regime (as galaxies approach and cross ~ L» in mass) in 
which today's galaxies transition from accretion-dominated, sec- 
ularly stable (dark matter-dominated) disks, to secularly unstable 
(self-gravitating) disks, which could quickly amplify ~ 10% ampli- 
tude perturbations into very strong bars and build significant pseu- 
dobulges, until later mergers destroy the remains of the disk and 
build massive classical bulges. 

This has important implications for the lifetimes of secu- 
lar processes of interest. The above comparisons assume disks 
where the internal response occurs over a single orbital period; 
if the systems have higher effective stability (i.e. secular re- 
sponses build more slowly), then the regime where they can be 
considered isolated for this time shrinks. Large-scale modes that 
require more than a few disk periods to self-amplify at low- 
redshift, or more than just a single disk period at high red- 
shift (z > 2), should be considered cosmologically dynamical 
systems (Figure [6]l - the galaxy grows comparably over this 
self-amplification timescale. Indeed, various observations of disk 
sizes and structure suggest that disks are sufficiently thick or 



are in this interesting range (Bell &de Jong 


1200 ll: iMcGaughl 


20051: ICourteau et alj 20071; iGilmore et alj|2002l; 


IWvse et alj|2006t 


Barteldrees & Dettmar 1994; de Griis et al.lll997 





Even if modes can evolve/self-amplify quickly such that a bar 
will grow efficiently and saturate at some final amplitude, these 
competing timescales define a maximum "isolated" lifetime for that 
saturated mode that is of interest, ~ O.lfHubbie (Figure [7J. There 
has been substantial debate regarding the lifetime of stellar bars 
in disks; but if modes live stably in isolation for longer than this 
time, they will encounter significant cosmological effects including 
e.g. significant new disk growth and mergers. Indeed, most studies 
do agree that lifetimes in isolation are at leas t this long (see e.g . 
Weinberdl 19851: iHernauist & Weinber j|l992bl : lFriedli et al.lll994l : 
Athanassoulall2002al : iKaufmann et alj .2007). Evolution of modes 
on longer timescales (e.g. some self-damping or buckling pro- 
cesses) should ideally be considered in a live cosmological con- 
text - the time in isolation may strengthen modes against exter- 
nal effects, but various studies have found that a moderate level 
of new gas accretion or passages of new substructure can dramati- 
cally change mod e evolution, both exciting and destroying bars and 
spiral waves (see Ath anassoula et al. 2005; Bournaud & Combed 



l2002l;lBerentzen et al JI2003I. |2004 120071 : iFovle et al.ll2008h ; not to 
mention that the presence of pre-existing strong bars may in turn 
affect these accretion/merger processes. 

Moreover, if modes live this long, their duty cycles will de- 
couple from those of their drivers. Even if, for example, all large- 
scale bars were initially driven by encounters with satellite galax- 
ies (minor interactions), if the isolated lifetime were much longer 
than this value, there would be no surviving correlation between the 
presence of bars and such companions. There has been consider- 
able observational debate regarding whether or not strongly-barred 
galaxies exhibit any strong preference for minor companions; cer- 
tainly th ere are at least many such galaxies w i thout close neigh- 
bors (seelElmegreen et alJl990l:IOdewahrJl994l : lMoles et aljl995l ; 
iMarquez & Moled ll996ULi et all l20oiT and references therein). 
This may in fact be because strong bars are not driven; however, 
it could also be consistent with the hypothesis that all such bars 
were initially driven, but are sufficiently long-lived. Constraints on 
bar lifetimes are needed to break the degeneracies. 

The level of cosmological dynamics also has implications for 
the numerical considerations involved in simulations of "isolated" 
systems. Properly following resonant self-interactions of bars may 



im ply steep resolution requ ire ments in N-body experiments (see 
e.g. Weinberg & Katzl2007al lb1: ICeverino & Klvpiij2007l : ISelIwoodl 
2008). However, there are other properties for which increasing the 
resolution in idealized cases may not be a more accurate represen- 
tation of reality. In terms of shot noise in the potential, for exam- 
ple, a model MW-like disk with 3> 10 6 particles will have potential 
fluctuations from smooth axisymmetry 5(j>/(j) < 1% over the spa- 
tial/timescales of interest (disk size and dynamical time). In cosmo- 
logical simulations, although the central regions of halos are rela- 
tively smooth, even dark-matter only simulations yield comparable 
or larger variations in the local po tential/velocity dis persion at e.g. 
MW-like disk effective radii (see IZemp et alj|2008f) . Even where 
smooth in space, such systems are not constant in time (as in ide- 
alized cases) at this level over several dynamical times. Moreover 
inclusion of baryons (which are not stripped efficiently, unlike dark 
matter subhalos which are efficiently destroyed at small radii and 
so do not "survive" to contribute substructure inside the centers 
of halos) enhances the clumpy, minor spatial substructure. In the 
Milky Way, for example, the LMC-SMC system represents a real 
deviation from a smooth, axisymmetric potential at a level larger 
than this limit near the solar radius. Ideally, tracking the evolution 
of substructure at higher resolution should involve not just a larger 
number of particles, but cosmologically motivated descriptions of 
substructure and accretion. 



Interestingly, at all redshifts, we find that traditional iso- 
lated "secular" evolution is most applicable around perturba- 
tions of fractional amplitude ~ 10%. This is a very interesting 
regime of parameter space: to the extent that it represents a frac- 
tional amplitude of substructure/accretion flows, it is a channel 
by which halos and low-mass galaxies gain much of their mass 
(e.g. Governato et alj|2007t iKazantzidis et a l. 2008; Stewart et al.1 
2008b). Moreover, "pseudobulges," associated with bulge forma- 
tion from secular evolution (e.g. bar-induced inflows and bar 
buckling; see e.g. lO'Neill & D ubinski 2 0ol; iMaver & Wadslevl 
120041 : iDebattista et al. 2004; Athanassoula 2005, and references 
therein) appear to dominate the bulge pop ulation at mass ratios 
of similar amplitude (B/T < 0.1 — 0.2; see Kuijken & Me rrifieldl 



1995l;|jogee et al.l2004l;lKormendv & Kennicutj20 04; Fishe ifeOOd : 
Fisher & Drory||2008l : IWeinzirl et alj|2009h . Suggestively, this also 
corresponds to typ ical amplitudes of observed strong bars (ref- 
erences above a nd Eskridge et al. 2 000l; Laurikainen et al. [J2002I : 



ISheth et alj|2003l ; iMarinova & Jogeell2007l ; lBarazza et alj|2 009) 



Of course, real systems exhibit more complex behavior then 
the simple scalings we derive here. Ultimately, detailed progress 
in modeling the interplay between continuous accretion of new 
substructure and cosmological driving of perturbations coupled to 
non-linear modes in galactic disks will require high-resolution N- 
body and hydrodynamic cosmological simulations. Some progress 
has begun towards mo deling these processes in a proper cosmo- 
logical context (see e.g. [Bournaud & Combe s 2002 ; Gauthie r"et al.l 



20061: iBerentzen & Shlosmarj 120061: IKaufmann et al.1 12001; 



Governato et al.l 2007; Foyle et al. 2008; Kazantzidis et al.1 



2008; 



Romano-Diaz et al.1 120081) - these studies highlight a key point 



here, that in a large regime of parameter space it is difficult to 
disentangle "secular" and cosmological processes. Our goal here is 
not to derive a rigorous quantitative description of one or the other. 
However, the simple arguments here should help to constrain and 
focus the discussion of where and when (in realistic cosmological 
settings) "isolated" evolution is important. 
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